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ABSTRACT 

The  application  of  phase  plane  analysis  to  a  missile  jet 
relay  control  system  is  reviewed,  explained,  and  illustrated. 
Only  a  single  lead  network  is  considered  under  conditions 
in  which  aerodynamic  damping  is  negligible.  The  assump¬ 
tion  that  all  forces  occur  as  step  inputs  is  used  to  simplify 
analysis  techniques  and  does  not  constitute  a  significant  limi¬ 
tation  on  the  application  of  this  method  for  design  purposes. 

To  aid  in  system  design,  an  analytic  method  for  the 
solution  of  phase  plane  plots  is  developed.  The  analytic 
technique  reduces  the  time  necessary  to  evaluate  a  system 
design  by  minimizing  the  necessity  for  construction  of  phase 
plane  diagrams. 

In  conclusion  a  summary  of  general  principles  which 
may  be  used  to  optimize  a  missile  jet  relay  control  system 
is  presented. 


PROBLEM  STATUS 

This  is  an  interim  report  on  one  phase  of  the  problem; 
work  is  continuing. 


AUTHORIZATION 
NRL  Problem  A02-18 
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INTRODUCTION 


The  jet  relay  control  system  has  wide  application  in  the  field  of  missile  control  and 
guidance.  A  jet  relay  control  system  is  generally  composed  of  a  gyro  reference  unit,  a 
control  amplifier,  pilot  and  slave  relays,  and  solenoids  which  actuate  reaction  jets.  The 
operation  of  this  system  may  be  described  by  considering  control  in  a  single  axis.  The 
gyro  output  signal,  which  indicates  missile  angular  deviation  from  a  reference  (gyro 
error),  is  applied  to  the  control  amplifier.  After  modification  by  the  amplifier  the  signal 
is  applied  to  a  balanced- output  amplifier.  The  output  of  this  amplifier  is  biased  so  that 
the  magnitude  of  the  error  signal  must  exceed  a  predetermined  value,  plus  or  minus,  to 
cause  pilot  relay  operation.  A  heavy-duty  relay  is  slaved  to  the  pilot  relay  and  when  the 
error  exceeds  the  amplifier  "dead  space,"  the  action  of  the  slave  relay  energizes  a 
solenoid  which  turns  on  a  reaction  jet.  The  action  of  the  jet  creates  a  correcting  moment 
which  rotates  the  vehicle  in  such  a  direction  as  to  reduce  gyro  error  and  cause  relay 
drop-out  (Fig.  1). 

This  type  of  control  system  is  most  effective  and  economical  where  relatively 
small  correcting  forces  are  necessary— usually  under  flight  conditions  wherein 
aerodynamic  disturbances  are  small  or  negligible. 

Since  this  method  of  control  is 
nonlinear,  frequency  response  anal¬ 
ysis  techniques  are  not  applicable, 
and  a  system  using  phase  plane  anal¬ 
ysis  has  been  generally  adopted. 

This  report  will  describe  a  phase 
plane  method  for  analysis  of  mis¬ 
sile  jet  relay  control  systems  and 
will  develop  analytic  design  criteria 
based  on  this  analysis.  The  tech¬ 
nique  described  in  this  report  makes 
the  simplifying  assumption,  valid 
for  design  purposes,  that  the  forces 
acting  on  the  vehicle  occur  as  step 
functions.  The  phase  plane  diagram 
is  a  plot  of  the  velocity  of  the  con¬ 
trolled  variable  as  a  function  of  its 
displacement.  This  technique  is  only 
valid  when  the  differential  equation 
describing  the  system  dynamics  may 
be  reduced  to  one  of  the  second  order. 

Missile  dynamics  in  the  absence  of 
aerodynamic  damping  meet  this 
criterion. 

An  example  of  a  phase  plane 
diagram  is  given  in  Fig.  2.  Con¬ 
sider  that  some  disturbance  of  short 
duration  has  resulted  in  an  angular  Figure  1  -  Typical  jet  relay  control  system 
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Figure  2  -  Simple  phase  plane 


velocity  80  of  the  missile.  The  angular  deviation  of  the  missile  from  the  zero  reference 
will  increase  with  time  and  will  generate  on  the  phase  plane  a  trajectory,  the  line  0,  1  and 
2.  When  the  angular  deviation  of  the  missile  has  reached  the  end  of  the  dead  space,  rep¬ 
resented  by  the  intersection  of  the  trajectory  and  the  pull-in  line,  point  1,  the  control 
relay  will  be  energized.  The  relay  actuates  a  solenoid  and  jet  thrust  is  developed.  Since 
there  is  a  time  delay  between  relay  energization  and  full  thrust,  the  angular  deviation  of 
the  missile  will  proceed  along  the  constant-velocity  path  during  this  time,  from  point  1  to 
point  2.  It  will  be  shown  later  that  when  the  reaction  jet  exerts  thrust  the  vehicle  trajec¬ 
tory  on  the  phase  plane  will  be  parabolic.  The  acceleration  due  to  the  jet  action  will 
reverse  the  velocity  along  the  parabola  from  point  2  to  point  3.  Since  relays  exhibit 
hysteresis,  the  angular  deviation  at  which  the  relay  is  de-energized  will  differ  from  the 
pull-in  line  as  represented  by  the  drop-out  line. 

There  will  be  a  time  delay  between  the  de- energization  of  the  relay  and  the  decay  of 
the  thrust,  during  which  the  trajectory  will  continue  along  the  parabola  from  point  3  to 
point  4.  The  system  then  has  obtained  a  constant  angular  velocity  opposite  in  sense  to 
the  initial  condition  and  the  sequences  will  be  repeated.  The  "limit  cycle"  is  established 
when  the  phase  plane  pattern  is  exactly  repeated. 

In  the  above  simplified  example  the  pull-in  and  drop-out  lines  were  represented  by 
vertical  lines  which  indicate  that  the  control  amplifier  did  not  modify  the  error  signal. 

In  practice,  lead  networks  are  employed  in  the  control  amplifier  in  order  to  damp  out 
transient  disturbances  and  to  increase  the  period  of  hunting  for  purposes  of  fuel 
conservation. 

This  report  will  develop  the  equations  of  the  relay  pull-in  and  drop-out  lines  for  a 
single  lead  network,  the  general  trajectory  of  vehicle  motion  on  a  phase  plane  plot,  the 
relationship  between  the  phase  plane  plot  and  time,  analytical  methods  for  determining 
the  system  limit  cycles,  and  general  rules  for  optimization  in  system  design. 


EQUATIONS  FOR  VEHICLE  TRAJECTORY  ON  THE  PHASE  PLANE 

From  an  examination  of  missile  dynamics  in  the  absence  of  aerodynamic  forces,  it 
can  be  seen  that  the  angular  deviation  8  of  the  missile  heading  with  respect  to  the  roll, 
pitch,  or  yaw  axis  is  given  by  the  following  equation: 

e  -  eo  +  eot  +  2  agt  ’  .  .  ,i  (1) 


and  differentiating  this  with  respect  to  t  gives 

8  ~  +  agt  ,  (2) 

where 

ag  =  the  system  angular  acceleration, 

80  =  the  initial  angular  displacement,  and 
80  =  the  initial  angular  velocity. 


By  combining  these  equations  to  eliminate  time,  the  angular  deviation  may  be  expressed 
as  a  function  of  angular  velocity: 


(g  -  gp); 

.  2a. 


(3) 


For  the  simplification  of  the  analysis,  all  quantities  should  be  nondimensionalized, 
since  then  a  parabolic  template  may  be  made  for  analyzing  the  system  for  various 
parameter  variations.  As  a  base  for  nondimensionalizing  all  quantities  having  time 
dimensions,  tb  =  1  second  is  used.  This  defines  a  nondimensional  time  scale  r  =  t/tb. 


To  nondimensionalize  the  other  quantities,  the  control  acceleration,  a,  is  used.  It 
follows  that  angular  displacement,  velocity,  and  acceleration  will  then  be  defined  as 


Hence  Eq.  (3)  may  be  written  in  the  nondimensional  form: 

^  =  ^0  +  a  tu  (^0  ~  V)  +  ~2  2  ^  ’  (4) 

which  may  be  simplified  to: 

n  =  no  ~  V)  •  (5) 

This  equation  represents  a  parabola  which  is  symmetric  about  the  fi  axis  and  with  vertex 
at 

n  =  °o  -  j  f  n02  . 
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There  are  four  general  conditions  which  may  govern  a  vehicle  trajectory  in  the  phase 
plane: 

(a)  No  disturbance,  no  correction,  ag  =  o.  In  this  case  the  vehicle  trajectory  is  rep¬ 
resented  by  a  constant-angular-velocity  line 

n  =  k. 

The  following  three  conditions  result  in  a  parabolic  trajectory: 

(b)  No  disturbance,  correction  applied,  ag  =  a.  The  parabola  will  be 

n  =  no  ±  |  -  ft02> . 

(c)  Disturbance,  no  correction,  ag  =  =  acceleration  due  to  the  disturbance.  The 

parabola  will  be 

«  =  no 

(d)  Disturbance,  correction  applied,  a  =  (aD+  a)  where  the  sign  of  a  will  be  opposite 
to  the  sign  of  fi  so  that  it  will  oppose  the  vehicle  angular  motion.  The  parabola  will  be 


Before  proceeding  with  examples  of  these  conditions,  however,  it  is  necessary  to  derive 
equations  for  the  relay  drop-out  and  pull-in  lines  and  develop  expressions  for  determina¬ 
tion  of  time  on  the  phase  plane. 


GENERAL  EQUATIONS  OF  RELAY  DROP-OUT  AND  PULL-IN  LINES 
(SINGLE  LEAD  NETWORK)  ON  THE  PHASE  PLANE 


To  determine  relay  drop-out  and  pull-in  lines  on  a  phase  plane  diagram  it  is  neces¬ 
sary  to  develop  an  expression  for  the  output  voltage  of  the  control  amplifier  as  a  function 
of  its  input  voltage.  This  expression  then  will  relate  the  controls  response  to  vehicle 
motion. 


For  most  relay  control  systems  the  control  amplifier  may  be  closely  approximated 
by  a  single  lead  network.  The  general  form  for  the  output  of  a  single  lead  network  in 
terms  of  LaPlace  transforms  is: 


where 


Eo(S)=ip.+TS)  E.(S) 

(1  +  rs) 


£„($■)  =  output  voltage, 


E.(S)  =  input  voltage,  and 
7,T  =  constants  of  lead  network. 

This  expression,  written  as  a  differential  equation,  becomes 


,dEe(t) 


dEi(t) 


(7) 


rEo(t)  +  T — Jj—  =  Ri(t)  +  T  dt 
This  is  a  first  order  linear  differential  equation  which  has  a  general  solution  as  follows: 


1  jj[  [*1(0  4eH  dtj  • 


E„(t)  =  e 

The  integration  may  be  performed  by  parts  and  the  equation  becomes 


(8) 


E0(t) 


JlE.(t)  +  —  IZHl}  F.ft'i  -  p~  T 

y  iv  '  y 


Ei(t) 


r 


t  y 

J  Ej(t)  e  T  *  dt  . 


(9) 


The  above  method  of  expansion  could  be  continued,  but  for  most  control  problems  a 
sufficient  degree  of  accuracy  may  be  obtained  by  assuming  that  all  vehicle  angular 
accelerations  occur  in  the  form  of  step  functions,  i.e.,  EA(t)  is  either  constant  or  zero. 
Thus  Eq,  (9)  becomes 


Since  the  input  voltage  is  directly  proportional  to  vehicle  angular  motion,  i.e., 
^i(t)  =  K6(t),  Eq.  (10)  maybe  rewritten  in  terms  of  9,  vehicle  angular  deviation: 


Ei(t)  +7  (7-1)  Ej(t)  -X-  (y-i)  E.(t) 


(10) 
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yE  (t)  >r  .  >r2  ..  /  ---pt 

— =  0r(O  =  0(0  +i(y- 1)  ^ft)  --—(7-1)  0(t)  ^1  —  e  T 


where  0R(t)  is  that  value  of  9  corresponding  to  the  relay  drop-out  or  pull-in  point  when 
there  is  no  lead  circuit  in  the  system.  Equation  (11)  is  the  exact  solution  for  the  relay 
pull-in  and  drop-out  lines  for  a  constant  angular  acceleration.  In  most  control  problems 
the  transient  term  may  be  considered  zero,  and  the  approximate  equation  for  the  pull-in 
and  drop-out  lines  may  be  used: 


0R(t)  =  B( t)  +T(y-  1)  0( t) 


(7-  1)  0(0  . 


Evaluation  of  the  transient  term  is  discussed  in  Appendix  A.  Since  the  transient  correc¬ 
tion  is  generally  insignificant  it  will  be  neglected  in  the  remainder  of  this  report. 

The  general  equation  for  relay  pull-in  and  drop-out  lines  may  be  nondimensionalized 
as  follows: 

„  „  .mfl  me  fi 

•  (I® 

where 


">  =  —  (7-  1)  . 


t.  =  1  sec  , 


Let  nRP  be  the  nondimensional  angular  position  for  relay  pull-in  determined  by  the  dead 
space.  Then  the  relay  drop-out  value  fiRD  is  given  by 


‘RD  -  “rp 


fi„p(R)  . 


where  K  is  the  value  of  the  relay  hysteresis. 

It  can  be  seen  that  the  positions  of  the  relay  pull-in  and  drop-out  lines  depend  on  the 
system  angular  acceleration;  this  important  point  will  be  demonstrated  in  an  example 
later.  Moreover,  since  the  system  angular  acceleration  is  considered  a  constant  at  any 
time,  the  relay  pull-in  and  drop-out  lines  are  represented  by  straight  lines  on  the  phase 
plane  plot,  with  slopes  m  and  intercepts  determined  by  fiR,  fi,  m  and  c. 

There  are  four  general  conditions  which  affect  the  equations  for  relay  pull-in  and 
drop-out: 

(a)  No  disturbance,  no  correction,  ay  =  0.  The  equation  for  the  relay  pull-in  line  in 
this  case  will  be 

±  fipp  =  fi  +  mfi  .  (15) 

(b)  No  disturbance,  correction  applied,  ap  =  a  (note  that  a  is  always  opposite  in  sign 
to  fiR).  The  equation  for  the  relay  drop-out  line  will  be 

±  fiRD  =  fi  +  mfi  ±  me  .  (16) 

(c)  Disturbance,  no  correction,  ag  =  an  (consider  >  o).  The  equation  for  the 
relay  pull-in  line  will  be 


±  fiRp  =  fi  +  mfi  +  me 


•m-rr-trH'V 
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(d)  Disturbance,  correction  applied,  aR  =  aD  +  a  (consider  aD  > 
the  relay  drop-out  line  will  be 

.  |txn  +  a| 

±  fiRD  =  n  +  mfl  ±  me  — — j - , 

where  the  sign  of  a  will  be  opposite  to  the  sign  of  nRP. 


I 


i 


o).  The  equation  for 

(18) 


RELATION  BETWEEN  THE  PHASE  PLANE  PLOT  AND  TIME 


Time  is  a  parameter  necessary  in  constructing  the  phase  plane  diagram,  but  it  does 
not  appear  explicitly  on  the  plot.  Therefore  a  relationship  between  time  and  the  phase 
plane  parameters  must  be  established. 

Again,  four  general  conditions  affect  the  relationship  between  time  and  the  phase 
plane  parameters. 

(a)  No  disturbance,  no  correction,  a  =  0.  Under  these  conditions,  the  trajectory  is 
represented  by  a  constant-angular- velocity  line,  Cl  =  K.  Then 


therefore 


I 


I 


At 


,  ah 
^ir 


(19) 


(b)  No  disturbance,  correction  applied,  ag  =  a.  In  this  case 

8  ^  =  8  j  +  a(t2~tj) 

=  8  j  +  a  A  t  , 


thus 


At=^. 


This  expression  is  then  nondimensionaiized  as  follows: 


“t. 


At  =  AO  tb  . 


(20) 


(c)  Disturbance,  no  correction,  a.?  =  aE.  In  this  case 

A 9  =  Oj,  At 


At  = 


A n  at. 


“D 


(d)  Disturbance,  correction  applied,  ag  =  +  a.  By  considerations  similar  to  the 

above, 


At  = 


AHtb  |  xx| 

I“d  +  «! 


(22) 
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EXAMPLE  PHASE  PLANE  ANALYSES 


The  use  of  the  equations  developed  for  phase  plane  analysis  can  best  be  demonstrated 
by  examples.  The  following  examples  are  not  intended  to  illustrate  design,  but  to  present 
technique. 


EXAMPLE  1.  SYSTEM  WITH  INITIAL  DISTURBANCE  ONLY 

Consider  the  case  where  a  control  system  is  subjected  to  an  initial  disturbance  which 
lasts  for  a  finite  time;  after  this  time  no  disturbance  is  present  (Fig.  3).  The  problem  is 
to  determine  the  maximum  excursion  experienced  by  the  system  due  to  the  disturbance, 
and  the  system  efficiency  in  the  final  limit  cycle. 

The  initial  disturbance  is  expressed  in  terms  of  a  disturbing  acceleration.  The  sys¬ 
tem  is  considered  at  rest  at  n  =  h  =  0  ;  and  n^,  nRD,  a,  m,  c,  and  the  system  time  delays 
are  known  quantities. 

Step  1.  Plot  on  the  phase  plane  the  pull-in  and  drop-out  lines. 

(a)  Since  the  initial  pull-in  may  occur  with  a  disturbance  and  no  correction, 

(Line  C)  nRP  =.n  +  mO  -  mc(^)  .  (23) 

(b)  Until  the  plot  has  been  started  it  will  not  be  known  whether  or  not  the  disturbance 
will  still  exist  at  drop-out.  Plot  the  drop-out  lines  for  no  disturbance. 


Figure  3  -  Phase  plane  with  initial  disturbance  only 
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(Lines  A) 


±  nRD  =  fi  +  mfi  ±  me 


(c)  Plot  the  pull-in  lines  needed  for  the  limit  cycle,  the  case  where  no  disturbance 
exists. 


(Line  B) 


±  nRP  =  fi  +  m  fi . 


Step  2.  Determine  the  length  of  time  the  initial  disturbance  will  exist  on  the  phase 
plane  and  plot  the  initial  trajectory. 

(a)  The  initial  trajectory  will  be  the  parabola 

o  -  l  a  62 

:  •  (26) 

(b)  This  parabola  will  exist  for  the  duration  of  the  disturbance  if  the  pull-in  line 
(line  C)  is  not  intersected  during  that  time.  The  time  duration  AtB  of  the  disturbance  is 
known;  the  parabola  will  extend  to  a  velocity  fix  defined  by  the  equation 

AtD  =  ^<V  (27) 

From  Fig.  3  it  may  be  seen  that  the  pull-in  line,  Line  C,  is  not  intersected  during:  the 
time  of  the  disturbance,  point  0  to  point  1.  When  the  disturbance  stops  the  pull-in  line 
shifts  and  becomes  line  B  since  there  is  no  acceleration.  In  this  case  pull-in  occurs  at 
the  instant  the  disturbance  stops,  since  the  trajectory  already  has  intersected  line  B;  he., 
pull-in  is  not  initiated  from  line  B  but  at  the  end  of  the  disturbance. 

Step  3.  Because  of  the  time  delay  in  relay  action,  solenoid  action,  and  jet  thrust 
buildup,  the  control  acceleration  will  not  occur  for  some  finite  time  after  the  initiation 
of  pull-in,  when  the  disturbance  stops.  This  time  delay  At,  is  known  for  the  system; 
during  this  time  delay  vehicle  angular  motion  is  along  a  constant-velocity  path,  point  1  to 
point  2.  The  length  of  this  path  is  determined  by  the  equation 

n2  -  n1  =  An  =  nxAt .  (28) 

Step  4.  At  point  2,  the  vehicle  angular  motion  is  due  to  the  control  angular  accelera¬ 
tion  a,  and  the  trajectory  in  the  phase  plane  may  be  expressed  as 

n  =  n2  -  \  (h2  -  n2) .  (29) 

In  practice  this  portion  of  the  trajectory  is  drawn  using  a  parabolic  template 

fi  =  j  n2 .  (30) 

This  template  is  moved  with  vertex  on  the  n  axis  until  the  curve  intersects  the  trajectory 
at  point  2.  The  trajectory  is  then  drawn  along  the  constant-acceleration  path,  point  2  to 
point  3. 

Step  5.  The  maximum  vehicle  angular  excursion  is  the  intersection  of  the  parabola, 
from  point  2  to  point  3,  with  the  fi  axis: 

(31> 


since  tb  =  1  second. 
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Step  6.  Point  3  is  the  intersection  of  the  trajectory  and  the  drop-out  line.  Again, 
during  the  time  delay  in  the  relay,  solenoid,  and  thrust  decay,  the  trajectory  in  the  phase 
plane  will  continue  along  the  parabolic  path.  The  point  at  which  thrust  is  zero,  point  4, 
is  determined  from  the  equation 

At  =  n4  -  n3  ,  (32) 

where  At  is  the  time  delay,  and  since  tb  =  1  second.  This  procedure  is  then  continued 
until  the  pattern  repeats  itself,  at  which  time  the  limit  cycle  has  been  reached. 

With  the  limit  cycle  a,  b,  c,  and  d  established,  it  is  possible  to  determine  the  system 
efficiency.  A  convenient  way  of  expressing  the  efficiency  is  by  the  ratio,  jet-on  time  per 
cycle  divided  by  the  period  of  cycle.  Hereafter  this  ratio  will  be  used  as  the  definition  of 
efficiency.  It  should  be  carefully  noted  that  the  higher  this  ratio,  called  the  "Efficiency 
Number,"  the  lower  the  efficiency  of  the  system: 


eff.  no. 


jet-on  time  per  cycle 
period  of  cycle 


(33) 


With  the  efficiency  established,  the  total  jet-on  time  for  any  flight  duration  may  be  deter¬ 
mined  by  multiplying  the  efficiency  by  the  flight  time. 


Step  7.  Determine  the  jet -on  time  by  taking  the  value  of  dQ  and  do  at  a,  b,  c,  or  d.  In 
the  limit,  the  pattern  should  be  symmetric  and  the  absolute  value  of  these  quantities 
should  be  the  same  at  any  of  the  four  points.  In  this  case 


jet-on  time 


2(2 n0) , 


(34) 


jet-off  time  = 


(35) 


EXAMPLE  2.  SYSTEM  WITH  A  CONTINUOUS  DISTURBANCE 

Consider  the  case  in  which  a  system  is  subjected  to  a  continuous  disturbance.  If  the 
acceleration  due  to  the  disturbance  is  large  enough,  the  trajectory  of  the  vehicle  in  the 
phase  plane  will  oscillate  about  the  positive  pull-in  and  drop-out  lines,  as  shown  in  Fig.  4. 
There  is,  however,  a  possibility  that  the  vehicle  trajectory  will  intersect  the  negative 
pull-in  line,  resulting  in  a  phase  plane  unlike  that  of  Fig.  4.  Because  of  this  possibility 
the  negative  pull-in  and  drop-out  lines  should  be  plotted.  The  problem,  as  in  Example  1, 
is  to  determine  the  maximum  excursion  experienced  by  the  system  due  to  the  disturbance, 
and  the  system  efficiency  in  the  final  limit  cycle. 

Again  the  disturbance  is  expressed  in  terms  of  a  disturbing  acceleration.  The  sys¬ 
tem  is  considered  at  rest  at  d  =  d  =  0;  and  dpP,  a,  m,  c,  and  the  system  time  delays 
are  known  quantities.  The  phase  plane  for  this  condition  is  plotted  as  follows: 

Step  1.  Plot  the  drop-out  and  pull-in  lines: 

dp  =  H  +  md  ±  m  c  d .  (36) 

(a)  For  the  pull-in  line  A,  d  is  dependent  on  the  acceleration  of  the  disturbance: 


ft  =  _£ 

pull-in  a 


(37) 
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Therefore 

±  ClR p  =  0  +  mfl  -  me  •  (38) 

Both  pull-in  lines  should  be  plotted  to  determine  the  possibility  of  the  vehicle  trajectory 
intersecting  the  negative  pull-in  line. 

(b)  For  the  drop-out  line  B,  fi  depends  on  the  acceleration  due  to  the  disturbance  plus 
the  correcting  moment  acceleration: 


where  aD  >  0  and  the  sign  of 
sign  of  nRr.  Therefore 


a. 


drop-out  a  * 

a  is  such  as  to  oppose  vehicle  motion,  i.e.,  opposite  to  the 


(39) 


I  aD  +  a 

±  =  fi  +  mfi  +  me  | - - -  . 

Step  2.  Determine  the  initial  trajectory  due  to  the  disturbing  acceleration, 
(a)  The  initial  trajectory  will  be  the  parabola 


(40) 


n  = 


2  ar 


(41) 


(b)  This  trajectory  will  continue  after  intersection  with  the  pull-in  line,  point  1, 
until  point  2.  This  length  is  determined  by 


An  =  At 


a  ’ 


(42) 


where  At  is  the  time  delay  of  the  system  resulting  from  relay  action,  solenoid  action,  and 
thrust  build-up. 
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Step  3.  At  point  2  the  vehicle  motion  wiil  be  the  result  of  the  disturbance  and  an 
opposing  control  acceleration: 


n2  =  fii  + 1  ^  &  -  n?) .  (43) 

This  trajectory  will  continue  after  intersection  with  the  drop-out  line  at,  point  3  until 
point  4,  because  of  the  time  delay  of  the  system: 


An  =  At 


(44) 


Step  4.  At  point  4  the  vehicle  motion  is  again  the  result  of  the  disturbance.  The  cycle 
is  now  repeated,  using  the  same  parabolic  equations,  until  a  limit  cycle  a,  b  is  reached. 
The  maximum  excursion  nmwill  occur  at  the  extreme  intersection  of  the  1>  axis  and  the 
parabola  representing  the  disturbance  plus  the  correction.  The  efficiency  number  is  the 
same  as  that  defined  in  Example  1: 


eff.  no. 


jet-on  time  per  cycle 
period  of  cycle 


(45) 


In  the  present  case 


jet-on  time  =  2fiUBlit  cycle|_£_|; 


jet-off  time-=  2nUmit  cycle|jq- 


(46) 

(47) 


1 


i 


E 


"M'f 


ANALYTIC  SOLUTIONS  FOR  PHASE  PLANE  PROBLEMS 


It  is  possible  to  determine  most  of  the  necessary  design  criteria  for  a  jet- relay  con¬ 
trol  system  analytically.  An  analytical  approach  to  the  solution  of  phase  plane  diagrams 
makes  it  possible  to  use  digital  computer  facilities  and  reduces  the  time  required  to 
evaluate  system  design  variations. 

In  the  following  pages,  equations  which  may  be  used  to  solve  for  the  necessary  design 
parameters  of  a  jet  relay  control  system  will  be  developed  and  the  limitations  of  this  type 
of  solution  will  be  explained. 


SYSTEM  WITH  NO  DISTURBANCE 

The  phase  plane  diagram  for  a  system  with  no  disturbances  will  be  of  the  form  shown 
in  Fig.  5.  This  is  a  diagram  showing  the  system  limit  cycle.  At  point  0,  jet  control 
thrust  becomes  effective,  resulting  in  a  parabolic  vehicle  trajectory  on  the  phase  plane 
which  intersects  the  relay  drop-out  line  at  point  1.  From  1  to  2,  the  vehicle  continues  on 
the  parabolic  path  for  a  distance  corresponding  to  the  drop-out  time  delay.  From  point  2 
to  point  (-3)  the  vehicle  has  a  constant  velocity  since  no  acceleration  exists.  At  point  (-3) 
the  vehicle  trajectory  intersects  the  relay  pull-in  line.  Vehicle  motion  continues  along 
the  constant  velocity  line  from  point  (-3)  to  point  (-0),  a  distance  corresponding  to  the 
pull-in  time  delay.  From  point  (-0)  where  control  thrust  again  becomes  effective,  the 
vehicle  motion  will  continue  in  a  manner  similar  to  that  previously  explained  and  return 
to  point  0.  It  may  be  shown  analytically  that  the  phase  plane  diagram  of  the  limit  cycle 
for  a  system  with  no  disturbances  will  be  symmetric  about  the  ft  and  ft  axes. 


Figure  5  -  Phase  plane  with  no  disturbance 
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Alteration  of  the  pattern  of  vehicle  motion  on  the  phase  plane  may  be  accomplished 
by  varying  the  lead  circuit  parameters  and  thus  changing  the  pull-in  and  drop-out  line 
positions,  by  changing  the  drop-out  or  pull-in  time  delays,  or  by  changing  the  control 
system  dead  spot.  It  may  be  possible  by  varying  these  parameters  to  over-compensate 
the  system,  in  which  case  the  system  will  become  semicontinuous. 

The  case  of  an  over-compensated  system  is  illustrated  in  Fig.  6.  At  point  0,  jet 
control  thrust  becomes  effective  and  results  in  a  parabolic  vehicle  trajectory  on  the 
phase  plane  which  intersects  the  relay  drop-out  line  at  point  1.  The  distance  from  point  1 
to  point  3  on  the  parabola  corresponds  to  the  drop-out  time  delay.  In  this  case  the  vehi¬ 
cle  trajectory  intersects  the  relay  pull-in  line,  point  2,  before  completion  of  the  relay 
drop-out  time  delay.  Under  this  condition  the  pull-in  time  delay  begins  at  point  2  and 
extends  to  point  (-0).  In  other  words,  part  of  the  pull-in  time  delay  is  occurring  during 
the  drop-out  delay.  At  point  3  the  drop-out  time  delay  ends  and  the  vehicle  moves  with  a 
constant  velocity  from  point  3  to  point  (-0),  a  distance  corresponding  to  the  remainder  of 
the  pull-in  time  delay.  The  pattern  from  point  (-0)  to  point  0  is  symmetric  with  the 
previous  pattern.  In  the  development  of  an  analytic  solution  for  the  limit  cycle  of  a  sys¬ 
tem  with  no  disturbances  the  over-compensated  system  is  not  considered,  since  this  type 
of  system  is  inefficient. 

Since  the  phase  plane  pattern  for  a  limit  cycle  (Fig.  5)  is  symmetric  with  respect  to 
the  n  and  fi  axes,  it  is  possible  to  completely  describe  the  properties  of  the  system  by 
knowing  the  coordinates  of  only  one  of  the  points  where  a  slope  discontinuity  exists. 
Therefore,  the  problem  of  solving  analytically  for  the  properties  of  a  system  in  the  limit 
cycle,  reduces  to  the  problem  of  solving  for  the  coordinates  fi0,  no  at  point  0.  It  is  known 
from  the  symmetry  of  the  pattern  that 

no  =  «2  =  "n-o  =  ~n-2>  (48) 

60  =  -n2  =  -n.„  =  n.2.  (49) 

The  first  step  is  to  solve  for  the  intersection  of  the  relay  drop-out  line  and  the  parabolic 
trajectory.  The  equation  for  the  relay  drop-out  line  will  be 


A 


Figure  6  -  Over-compensated  system 
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fiRD  ~  flj  +  mHj  +  me  , 

and  the  equation  for  the  parabolic  trajectory  will  be 


<V=  no  -  "AV 


Solving  these  equations  simultaneously  yields 


r~ — ~ 

(  h2\ 

m2  -  2 

(0^-  me)  ~(no+-f-) 

=  m  ± 


where  the  minus  sign  applies  in  the  present  case.  Now, 

+  |  fij  -  flj|  =  AtD  , 

where  AtD  is  the  time  delay  at  drop-out;  since  by  Eq.  (49)  n2  =-  fi01 


n  -  Atn  —  n, , 


and 


H0  -  AtD  +  m  =  1/  m2  -  2 


-  ^n0  +% 

Squaring  both  sides  of  this  equation  and  combining  terms  yields 

AtD2 

no  =  -y-  +  nofm  -  AtD)  -  At^n  +  OrD  -  me  . 

The  time  delay  at  pull-in  is 

n  -  n. 


Atp  ~  0 


and,  at  the  intersection  of  the  constant-velocity  trajectory  and  the  pull-in  line, 

^3  “  ^RP  “  "A  * 

combining  these  two  equations  yields 

=  ^RP  +  ^o^tP  ~  m)  ‘ 

This  egression  for  no  may  be  equated  to  Eq.  (56)  to  yield  a  solution  for  no: 

2 

'  r  ,  At-f 

no[2m  -  (Atp  +  Atp)]  = - +  At^  +  (fiRp  -  fiRD)  +  me  . 

For  the  case  where  the  time  delays  are  equal,  i.e.,  AtD  =  Atp, 

At2 


and 


n02(m  -  At)  =  -  -y-  +  At  m  +  (0^  -  nRD)  +  me 
fi0  =  (HRp  -  m  ft0  +  Atpfl0)  . 


(50) 

(51) 

(52) 

(53) 

(54) 

(55) 

(56) 

(57) 

(58) 

(59) 

(60) 

(61) 

(62) 
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When  the  coordinates  fi0,  fiD  of  point  0  are  known,  the  maximum  excursion  in  the  limit 
cycle  may  be  determined  from  the  parabolic  relationship 


o  =n+4-  n2 


(63) 


m  “  o  2  ° 

The  efficiency  number  of  the  system  may  be  determined  from  fi0  and  fio  as  follows: 

jet-on  time  =  2(2  fi0)  ,  (64) 


jet-off  t ime  =  2 


n 


period  =  4  fiQ  +  -r2-  , 


e  f  f .  no .  =  - 


1  + 


(65) 

(66) 

(67) 


It  should  be  reiterated  that  the  above  equations  apply  to  a  control  system  with  no 
disturbances  present,  and  that  the  equations  are  not  valid  for  an  over-compensated  system. 
One  indication  from  the  analytic  solution  that  a  system  has  been  over -compensated  is  the 
appearance  of  a  negative  no .  It  is  suggested  that  whenever  the  analytic  solution  is  in 
doubt  a  quick  graphical  check  be  made.  A  quick  means  for  a  graphic  check  is  to  plot  fi0 
and  no  at  all  four  points,  connect  the  trajectories  and  then  measure  the  delay  times.  The 
measured  times  should  be  equivalent  to  the  known  times. 


SYSTEM  WITH  A  CONTINUOUS  DISTURBANCE 

The  phase  plane  diagram  for  a  system  with  a  continuous  disturbance  will  be  of  the 
form  shown  in  Figure  7.  This  is  a  diagram  showing  the  system  limit  cycle.  At  point  0 
the  control  thrust  becomes  effective  and  the  result  is  a  parabolic  trajectory  for  which  the 
resultant  acceleration  is  due  to  the  difference  between  the  disturbance  and  the  control 
force.  The  vehicle  trajectory  intersects  the  drop-out  line  at  point  1.  From  there  to 
point  2  the  vehicle  motion  continues  on  the  same  parabolic  path  for  a  distance  corre¬ 
sponding  to  the  drop-out  time  delay.  From  point  2  to  point  3,  the  vehicle  trajectory  in 
the  phase  plane  is  again  parabolic  with  an  acceleration  due  to  the  disturbing  force  only. 
At  point  3,  the  trajectory  intersects  the  relay  pull-in  line.  The  vehicle  motion  continues 
on  the  parabolic  trajectory  due  to  the  disturbance  from  point  3  to  point  0,  a  distance 
which  corresponds  to  the  relay  pull-in  time  delay.  The  pattern  is  then  repeated. 

In  the  case  of  a  continuous  disturbance  the  phase  plane  pattern  in  the  limit  cycle  is 
symmetric  with  respect  to  the  fi  axis.  Because  of  this  symmetry  it  is  possible  to 
describe  the  properties  of  the  system  completely  when  the  coordinates  at  the  intersec¬ 
tion  of  the  two  parabolas  are  known. 

Because  of  the  symmetry  of  the  pattern  (Fig.  7) 

no  =  -  n2 ,  (68) 

nD  =  n2 .  (69) 

The  first  step  is  to  solve  for  the  intersection  of  the  trajectory  with  no  correcting  force 
applied,  and  the  pull-in  line.  The  equation  for  the  pull-in  line  will  be 
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a 


Figure  7  -  Resultant  phase  plane  for  a  system  with 
a  continuous  disturbance 


n3  =  o2  +  (n32  -  n22) .  (71) 

Solving  these  equations  simultaneously  gives 


The  next  step  is  to  solve  for  the  intersection  of  the  trajectory  with  correcting  force 
applied  and  the  drop-out  line.  The  equation  for  the  relay  drop-out  line  will  be 


Solving  these  equations  simultaneously  gives 


Next  the  time  delays  are  considered.  The  drop-out  time  delay  will  be 
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Since  a  and  will  always  be  opposite  in  sign  the  equation  may  be  rewritten  without  the 
absolute  value  signs  as  follows: 


period  of  cycle  =  2no  — +  ^ 


20, 


_ V 

\aD^a  ~  aD>/’ 


20 


a 


20„ 


°  (a  -  aD)  _  aD 
a2 


°  UD(a  -  aD) 


(87) 


eff.  no. 


a 


(88) 


SUMMARY  OF  EQUATIONS  FOR  THE  ANALYTIC  SOLUTION  OF 
PHASE  PLANE  PROBLEMS 


SYSTEM  WITH  NO  DISTURBANCE 

For  a  system  with  no  disturbance  the  characteristic  coordinates,  derived  originally 
as  Eqs.  (59)  and  (60),  are 


At  tT 

~  (^D  *  -  2~  *  ^  ( Opp  ~  Urd)  *  mc 

fiQ  =  Ojjp  +  no(Atp  —  m)  ; 

the  efficiency  number  of  the  system,  derived  originally  as  Eq.  (67),  is 

1 


eff.  no.  =■ 


1  4 


(89) 

(90) 


(91) 


and  the  maximum  excursion  in  the  limit  cycle  is 


n  =  n  4  i  rt 2  . 
m  o  2  ° 


(92) 


SYSTEM  WITH  A  CONTINUOUS  DISTURBANCE 

For  a  system  with  a  continuous  disturbance  the  characteristic  coordinates  are 


0o  [2m  ~  (AtD  4  Atp)] 


Atp  ja  4  aD[  Atp  |  Ogl 

2  |aj  2  lal 

|a  4  aD|  I Ojj I 


+  m  At 


(93) 


+  AtP  a  j  +  °RP  ~  ^RD  +  mc  ’ 


.  aD  _ 

a 

4  mAtp  —  +  fipp  +  mc 

a 

Atr 


the  efficiency  number  of  the  system  is 


eff.  no.  =  —  ; 

and  the  maximum  excursion  in  the  limit  cycle  is 


no(m  -  Atp)  ; 


(94) 


(95) 


n  =  nQ  4  .  |a| — r  Ci 2 

m  °  2|a  4  opl  ° 


(96) 
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GENERAL  PRINCIPLES  FOR  OPTIMIZATION  IN  SYSTEM  DESIGN 


The  design  of  a  relay  jet  control  system  consists  of  determination  of  the  control 
thrust  and  amplifier  dead  spot  necessary  to  satisfy  the  system  limitations  and  determina¬ 
tion  of  the  lead  network  parameters  which  will  result  in  the  greatest  system  efficiency 
(the  lowest  efficiency  number).  Since  every  system  will  have  some  limiting  angular 
excursion,  the  choice  of  the  amplifier  dead  spot  and  control  acceleration  will  be  based  on 
this  limitation. 


Consider  a  very  simple  system  with  a  disturbance,  no  lead  network,  no  time  delays, 
and  no  relay  hysteresis  (Fig.  8).  In  this  case  the  vehicle  trajectory  will  immediately  be 
in  a  limit  cycle  and  by  solution  of  the  equations  for  the  trajectory,  the  system  maximum 
excursion  may  be  expressed  as  follows: 


where 


n 


ip 


aD  \ 

O  -  aD) j  ’ 


(97) 


nD  =  amplifier  dead  spot, 
aD  =  disturbing  acceleration, 
a  =  control  acceleration,  and 


n  =  maximum  excursion. 

m 


The  consideration  of  time  delays  and  relay  hysteresis  would  increase  the  maximum 
excursion,  while  the  use  of  a  lead  network  would  reduce  it.  Since  these  effects  tend  to 


a 


Figure  8  -  Phase  plane  without  lead  circuit  or 
hysteresis 
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compensate  they  are  neglected  in  the 
initial  design  choice  of  control  accel¬ 
eration  and  dead  spot.  The  relation¬ 
ship  between  maximum  excursion,  con¬ 
trol  acceleration,  and  dead  spot  for  the 
simple  system  described  above  is 
shown  'in  Fig.  9. 

For  a  system  with  an  initial  dis¬ 
turbance,  which  eventually  oscillates 
in  a  limit  cycle  with  no  disturbance, 
the  efficiency  number  as  defined  earlier 
becomes  less  (the  efficiency  becomes 
greater)  with  increasing  a  and  also 
with  larger  dead  spots. 

For  a  system  with  a  continuous 
disturbance  the  most  efficient  system 
has  the  largest  possible  control  accel- 


Figure  9  -  Dead  spot  vs  ratio  of  disturbing 
to  control  acceleration  for  different  maxi¬ 
mum  excursions 


eration,  since  in  this  case  efficiency  depends  only  on  control  acceleration.  The  optimum 
lead  network  will  depend  on  the  time  delay  and  the  relay  hysteresis  as  well  as  the  control 
acceleration  and  dead  spot. 


After  selection  of  a  control  acceleration  and  amplifier  dead  spot  which  satisfy  the 
system  limitations,  it  is  necessary  to  determine  the  optimum  lead  network. 

For  a  system  with  only  an  initial  disturbance,  the  lead  network  is  optimized  on  the 
basis  of  maximum  system  efficiency  in  the  limit  cycle.  Figure  10  shows  the  relationship 
between  7  and  T,  the  lead  circuit  parameters,  and  efficiency  number  in  the  limit.  From 
consideration  of  this  figure  the  value  of  y  chosen  for  the  lead  network  should  be  the 
largest  value  which  can  be  attained  and  still  provide  an  acceptable  system  noise  level. 
After  the  determination  of  y,  a  series  of  calculations  employing  the  equations  developed 
earlier,  may  be  used  to  determine  the  optimum 
value  of  T . 

Fora  system  with  a  continuous  disturbance, 
the  lead  network  is  optimized  on  the  basis  of  the 
smallest  maximum  excursion  in  the  limit  cycle. 

The  relationship  between  y  and  T,  and  the  max¬ 
imum  excursion  in  the  limit  cycle  are  shown  in 
Fig.  11.  After  the  largest  y  consistent  with 
system  noise  considerations  has  been  chosen, 

T  may  be  determined  by  a  series  of  calculations 
employing  the  equations  developed  earlier.  An 
interesting  side  light  on  this  type  of  system  is 
illustrated  in  Fig.  12.  This  figure  compares 
the  initial  system  excursion  (the  first  excursion 
experienced  by  the  system  if  it  is  initially  at 
rest  and  the  disturbance  appears  as  a  step  input) 
with  the  maximum  excursion  in  the  limit.  It 
may  be  seen  that  the  optimum  value  of  T  on  the 
basis  of  the  maximum  limit  cycle  excursion 
results  in  an  initial  excursion  which  has  approx¬ 
imately  the  same  value,  indicating  that  the  sys¬ 
tem  is  nearly  critically  damped  for  this  value 
of  T. 


EFFICIENCY 


Figure  10  -  Lead  circuit  parameter 
Tvs  efficiency  number  for  different 
values  of  y 
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■5  +  20 

-8  +  30  i  7+10  y 


Figure  11  -  Effect  of  different  values  of  y  on 
the  relation  between  the  lead  circuit  parame¬ 
ter  T  and  the  maximum  excursion  in  the  limit 
cycle,  for  a  system  with  a  constant  disturbance 


The  system  time  delay  has  not 
yet  been  mentioned.  Figure  13 
shows  the  relationship  between  the 
system  efficiency,  or  maximum 
excursion  in  the  limit,  and  the  sys¬ 
tem  time  delay.  It  is  obvious  that 
time  delays  should  be  minimized  as 
much  as  possible. 

Relay  hysteresis  generally  has  a 
relatively  small  effect  on  the  system 
efficiency. 
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then 

-5-  ind  -  y)  =  y  -  x.  (A7) 

For  any  value  of  the  lead  circuit  parameter  c/m  =  (l/y-  l),  x  will  approach  infinity 
as  y  approaches  unity.  Hence,  if  a  curve  of  y-x  =  Jn(i  -  y),  in  other  words  c/m  =  1  or 
y  -  2,  is  plotted  in  the  xy  plane  (Fig.  Al)  there  will  be  a  point  somewhere  on  this  curve 
with  the  same  value  of  y  as  would  exist  on  a  curve  where  c/m  has  a  value  other  than  unity. 
The  value  xD  =  fi0  -  fipp/mn  is  known  in  each  particular  case. 


X, 


Figure  Al  -  Comparison  of  transient  correction 
curves  (no  acceleration)  for  c/m  =  1 


It  is  required  to  find  yj  for  the  case  when  c/m  is  not  unity.  But  yt  will  also  be  on 
the  curve  of  y-  x  =  ln(l  -  y)  with  a  corresponding  abscissa  xlt  Substituting  the  coordi¬ 
nates  of  these  points  in  the  equations  of  the  two  curves  yields 


(A8) 

(A9) 

(A10) 

(All) 

(A  12) 

» 

(A13) 

(A14) 

(A15) 


x-intercept  =  (y  -  l)x 
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This  straight  line  may  be  drawn,  since  xQ  and  y  are  known  in  each  particular  case.  The 
intersection  of  this  line  with  the  curve  y-x  =  ln(l  -  y)  then  gives  the  required  value  of  y. 


Since  the  slope  of  the  straight  line  only  depends  on  the  value  of  y,  for  a  particular 
lead  circuit  all  lines  corresponding  to  different  values  of  xo  will  be  parallel  to  each  other 
The  procedure  for  finding  the  actual  relay  pull-in  point  when  fip'p,  y,  and  xD  -  Clq  -  f^p/ml  Q| 
are  given  is  therefore  as  follows:  On  Fig.  A2,  draw  the  straight  line  of  Eq.  (A12).  Read 
off  the  value  of  y  where  the  straight  line  intersects  the  curve  y-x  =  ln(l  -  y).  The  actual 
relay  pull-in  point  is  then  given  by 


_  ~ 
mj  n| 

f!p  =  m|  fi|  y  +  fipp  . 


(A15) 

(A16) 


Figure  A2  -  Transient  correction  curve  (no  acceleration)  for  c/m  =  1 


Operations  in  the  upper  half  of  the  phase  plane,  where  n  is  positive,  can  be  analyzed 
in  a  similar  manner.  Here  it  can  be  shown  that  if 


x 


m|n| 


^RP  ~  ^P 

m|n| 


the  foregoing  procedure  for  finding  Hp  may  be  used  except  that  now 

^P  =  ^RP  ~ ^  m  |  fl  |  y  . 


(A17) 

(A18) 


(A19) 
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p 


t  » 


For  a  lead  network  where  7  =  2,  it  is  only  required  to  calculate  x0;  then  y  is  directly 
read  off  the  curve  y  -  x  =  ln(  1  -  y). 


Correction  when  a  Constant  Acceleration  Exists 

To  determine  the  actual  relay  drop-out  or  pull-in  point  when  an  acceleration  exists 
consider  the  general  equation  for  relay  pull-in  or  drop-out 

=  ft  +  mn  -  men  (1  -  e  c).  (A2°) 

The  steady- state  drop-out  line  will  be 

ftp  =  ft  +  mft  —  me  ft  .  (A21) 

Let  ftA  be  the  approximate  relay  drop-out  point  (steady-state),  and  ftj.  the  actual  relay 
drop-out  point;  then 


ftp  —  ftA  =  mftA  -  me  ft  , 

ftR  -  ftj  =  mftj.  —  me  ft  ^1  -  e  c  j  . 
Dividing  both  equations  through  by  me  ft  gives 

ftp  -  ftj.  fty  ^ 

-  1  . 


meft  eft 


meft  eft 


(A22) 

(A23) 

(A24) 

(A25) 


Now  t ,  the  nondimensional  time  to  go  from  the  beginning  of  the  acceleration  8  to  the 
drop-out  line  will  be 


therefore 


A  ft  a  _  A  ft 

0  ~T 


ftp  -  Oj,  ftj 

- - —  =  -=--(!  -  e 

meft  eft 


Aft 

cfo 


(A26) 


(A27) 


In  order  to  solve  the  equations,  consider  the  phase  plane  axis  to  have  been  transformed 
so  that  ftj  =  0.  Then 


Aft  =  n2 . 


With  this  transformation,  substitute 


ftc 

ft 

■rr—-~  • 

ftmc 


=  y  1 


(A28) 

(A29) 

(A30) 


1 
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Then  Eqs.  (A25)  and  (A27)  become 


where 


(C  +  1)  -  x  =  y , 

(C  +  1)  -  x  =  y  +  e'y  , 


(A31) 

(A32) 


me  ft 


(A3  3) 


The  general  equation  for  the  trajectory  on  the  phase  plane  with  an  acceleration  ft 
will  be 


1  U2 

a,  -  n,  =  — ~ — . 
2  1  2  n 


Dividing  this  equation  through  by  men  gives 


n2  1  c  ^2 


Now  substitute 


and  the  equation  becomes 


where 


mcO  ^  m  c2ft2 

n2 

x  =  — —  , 
meft 

ft2 

y  _  eft  ’ 


x  "  D  =  y2  , 
l  m 


(A3  4) 

(A3  5) 

(A3  6) 
(A3  7) 

(A3  8) 


D=-4.  (A3  9) 

meft 

If  these  curves  are  plotted  it  may  be  seen  that  actual  relay  pull-in  or  drop-out  occurs 
before  the  trajectory  intersects  the  steady-state  pull-in  or  drop-out  curve  (Fig.  A3). 

Examination  of  the  equation  for  the  trajectory,  Eq.  (A38),  indicates  that,  plotted  on 
log- log  paper,  it  is  the  equation  of  a  family  of  straight  lines  of  constant  slope.  Hence  if 
Eqs.  (A31)  and  (A32)  are  also  plotted  on  log- log  paper,  the  actual  relay  drop-out  or  pull- 
in  point  will  be  given  by  the  intersection  of  the  straight  line  (representing  the  trajectory) 
drawn  through  the  approximate  relay  drop-out  or  pull-in  point  determined  from  the 
phase  plane  plot,  and  the  line  representing  the  actual  relay  drop-out  or  pull-in  point. 

In  summary,  the  procedure  for  finding  the  actual  relay  pull-in  or  drop-out  point 
when  an  acceleration  exists  is  as  follows: 

1.  On  the  phase  plane  plot  draw  the  steady-state  relay  pull-in  and  drop-out  lines. 

2.  Find  the  nondimensional  change  in  velocity,  A  ft,  from  the  initial  point  on  the 
parabola  to  the  approximate  relay  pull-in  or  drop-out  point  as  given  by  the  intersection 
of  the  parabola  and  the  steady- state  relay  pull-in  or  drop-out  line. 
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Figure  A3  -  Relation  between  actual  relay  pull-in 
and  drop-out  and  phase  plane  trajectory 

3.  Divide  |Afi|  by  ch. 

4.  Enter  this  value  |An|/cfi  on  the  transient  correction  plot  (Fig.  A4)  and  find  the 
corresponding  point  on  the  steady- state  line. 

5.  Through  this  point,  draw  the  straight  line  representing  the  parabola.  This  is  a 
straight  line  parallel  to  the  straight  line  labelled  as  parabola  (reference)  on  the  plot. 

6.  The  intersection  of  this  line  with  the  actual  drop-out  line  determines  the  actual 
drop-out  point.  Read  off  this  value  of  lAn/ciil . 

7.  Multiply  this  value  by  cfi  and  adjust  the  value  of  |Afl|  on  the  phase  plane  plot  to 
obtain  the  actual  relay  pull-in  or  drop-out  point. 


Figure  A4  -  Transient  correction  for  relay  drop-out  or  pull- 
in  line  when  an  acceleration  exists 


Appendix  B 
LIST  OF  SYMBOLS 


a  =  angular  acceleration  due  to  a  control  moment 
otp  =  angular  acceleration  due  to  a  disturbance 
ag  =  system  angular  acceleration 
y  =  lead  circuit  parameter 
6  =  vehicle  angular  displacement 
6  =  vehicle  angular  velocity 
8  =  vehicle  angular  acceleration 
0.o  =  initial  angular  displacement 
eo  =  initial  angular  velocity 

8-p  =  value  of  6  corresponding  to  the  relay  drop-out  or  pull-in  point  when  there  is  no 
lead  circuit  in  the  system 

r  =  t/tb  =  nondimensional  time  scale 

fi  =  nondimensional  angular  displacement 

H  =  nondimensional  angular  velocity 

=  nondimensional  angular  acceleration 

flp  =  amplifier  dead  spot 

nm  =  maximum  excursion  in  limit  cycle 

nRP  =  nondimensional  pull-in  point  determined  by  dead  spot 

nRD  =  nondimensional  drop-out  point  determined  by  dead  spot 

c  =  T/y 

E0(s)  =  output  voltage 
E.fs)  =  input  voltage 
m  =  T(y  —  l)/y 
t  =  time 

tb  =  1  second  (time) 

At  =  system  time  delay 
AtD  =  system  time  delay  at  drop-out 
Atp  =  system  time  delay  at  pull-in 
T  =  lead  circuit  parameter 

*  *  * 
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